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Abstract 

We consider high-energy quasifree single- and two-proton knockout reactions induced by electrons 
and protons and address the question what target-nucleus densities can be effectively probed after 
correcting for nuclear attenuation (initial- and final-state interactions). Our calculations refer to 
ejected proton kinetic energies of 1.5 GeV, the reactions {e,e'p), {jiPp) and {p,2p) and a carbon 
target. It is shown that each of the three reactions is characterized by a distinctive sensitivity to the 
density of the target nucleus. The bulk of the (7,pp) strength stems from the high-density regions 
in the deep nuclear interior. Despite the strong attenuation, sizable densities can be probed by 
{p, 2p) provided that the energy resolution allows one to pick nucleons from s orbits. The effective 
mean densities that can be probed in high-energy (e, e'p) are of the order of 30-50% of the nuclear 
saturation density. 
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Nucleon knockout studies in quasifree kinematics belong to the most powerful instru- 
ments for studying the structure of nuclei. Since the 1960's the A{e, e'p) reaction has pro- 
vided a wealth of information about the merits and limitations of the nuclear shell-mo del m. 

n 

Quasifree proton scattering from nuclei A(p, 2p) has a somewhat longer history [21] and could 
in principle provide similar information as A{e,e'p). With three protons subject to attenu- 
ation effects, in A{p, 2p) the description of the initial and final-state interactions (ISIs and 
FSIs), is a more challenging issue than in A{e,e'p). Recent applications of A[p,2p) involve 
the analyzing power (Ay) as an instrument for probing medium modifications of hadron 
properties and the density dependence of the nucleon-nucleon interaction . In inverse 
kinematics (i.e. the p{A, 2p)A — 1 reaction) the (p, 2p) process offers great opportunities 
for systematic studies of the density and isospin dependency of single-particle properties in 
unstable nuclei at high-energy radioactive beam facilities j^. Experiments of that type 
have the potential to study the equation-of-state for nuclei far from equilibrium. With regard 
to quasifree A{e,e'p), recent developments include the search for medium modifications of 
electromagnetic form factors through double polarization experiments of the type ^He(e, e'p) 
[3]. Another line of current research is that of the two-nucleon removal reactions A{e,e'pp) 
and A{'j,pp) in selected kinematics. These reactions provide a window on the short-range 
structure of nuclei [Sj. 

The development of an appropriate reaction theory is essential for reliably extracting the 
physical information from the nucleon knockout reactions. For nucleon kinetic energies up to 
about 1 GeV the distorted wave impulse approximation (DWIA) has enjoyed many successes 
in that it could reproduce a large amount of measurements fairly well Constraining the 
optical potentials for DWIA calculations, however, heavily depends on the availability of 
elastic proton-nucleus scattering data. Moreover, the optical potentials exhibit a substantial 
kinetic-energy dependence. This energy dependence makes it difficult to make more general 
statements about e.g. the role of attenuation effects and the effective densities that can be 
probed in the various reactions. At sufficiently high nucleon energies the Glauber approach 
provides a valid and highly efficient alternative for the DWIA framework. The Glauber 
approach has the advantage that the effect of ISIs and FSIs can be computed from the 
knowledge of the elementary proton-proton and proton-neutron differential cross sections 
and of the density of the target (residual) nucleus. Moreover, for nucleon momenta exceed- 
ing about 1 GeV, the energy dependence of the parameters entering the Glauber calculations 



2 



is relatively smooth. This results, for example, in measured and computed nuclear A{e, e'p) 
transparencies which exhibit little energy dependence nucleon kinetic energies larger than 

n 

0.5 GeV [9]. From the theoretical point of view, it allows one to make more universal state- 
ments about the predicted role of nuclear attenuation. Another advantage of the Glauber 

reactions, including electromagnetic and 
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approach is that it is applicable to a wide range of 
hadronic probes, with stable and unstable nuclei 

We exploit the robustness of the Glauber approach to study the density dependence of 
quasifree nucleon removal reactions. Indeed, investigations into the medium dependence of 
nucleon properties and the study of the nuclear structure of unstable nuclei e.g., heavily 
rely on the possibility of effectively probing regions of sufficiently high density in the target 
nucleus. Nuclear attenuation effects on the impinging and ejected protons can cause the 
nucleon knockout reactions to effectively probe regions of relatively small density near the 
surface of the target nucleus. The description of nuclear attenuation brings in a certain 
degree of model dependence. We stress the importance of making cross checks over different 
types of reactions (electromagnetic versus hadronic probes) and linking single-nucleon to 
two-nucleon knockout reactions. Here, we report on a study of the effective nuclear density 
that can be probed in reactions that have one nucleon {A{e,e'p)), two nucleons {A{'~f,pp)) 
and three nucleons {A{p, 2p)) subject to nuclear attenuation effects. 

In Ref. jl^ a relativistic extension of the Glauber method was introduced. The method 
was coined RMSGA (relativistic multiple-scattering Glauber approximation). In line with 
the assumptions of a typical relativistic DWIA model, the RMSGA model uses a relativistic 
mean-field to describe the target and residual nucleus in combination with the impulse 
approximation for the interaction Hamiltonian. The RMSGA differs from the relativistic 
DWIA in that it uses a relativistic extension of the Glauber method to treat initial-state 
and final-state interactions. 

In a factorized approach, the differential cross sections for the single-nucleon removal 
reactions considered here (i.e. A{p,2p) and A{e,e'p)) are proportional to the distorted 
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FIG. 1: (color online) The function 5{r,9) for the ^^C(e,e'p) and ^^C(p, 2p) reaction. For both 
types of reactions we consider an energy transfer of 1.5 GeV and a three-momentum transfer q 
that is tuned to probe the maximum of the momentum distribution (i.e. Pm=0 MeV for knockout 
from the sl/2-orbit and Pm=115 MeV for removal from the p3/2-orbit). For the {e,e'p) results 
the proton is detected along the direction of the momentum transfer. For the {p, 2p) the incoming 
proton has a kinetic energy of about 3 GeV and the two ejected protons have a kinetic energy of 
1.5 GeV. They are detected under an angle of about 32° but on opposite sides of the incoming 
beam. For the sake of reference, the proton root-mean-square radius in ^^C as determined from 
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elastic electron scattering is (r^) = 2.464 it 0.012 fm 
momentum distributions pfnK)(Pm.) 
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dr j deS (r, 9) , (1) 

where the quantum numbers {nnm) determine the orbit of the struck nucleon, (puKmi'f^ is the 
corresponding relativistic single-particle wave function and u{p, s) a four-component free- 
particle spinor. The missing momentum pm is determined by the difference between the 



D{r,e)= / d(f)r'^ sme——r-u{prn,s)Sly^sQj^{f)(p^^ (2) 



asymptotic three-momentum of the ejected nucleon p and the three-momentum transfer q. 
We define the 2;- axis along the q and the xz-plane as the reaction plane. The function 6{r, 6) 
defined in Eq. (IT]) encodes the contribution from an infinitesimal interval around r and 6 to 
a single- nucleon removal cross section [3j. The function D{r, 9) which was introduced in ([T]) 
reads 

-iPm-f 

The Glauber phase operator »5p^msga(^) encodes the combined effect of the initial and final- 
state interactions. Retaining only central interactions, which is a fair approximation at 
the kinetic energies considered here, it can be considered as a scalar operator. For the 
A{p, 2p) reaction, the (Srmsga becomes a mult i- dimensional convolution over the squared 
wavefunctions of the spectator nucleons times the profile functions for the impinging proton 
and two ejected protons [l4|. For the A{e,e'p) case the convolution involves the squared 
wavefunctions of the spectator nucleons times the profile function for the ejected nucleon 



We now wish to formulate the analog of Eq. ([T]) for two-nucleon knockout reactions of 
the type A{e, e'pp) and A{'y,pp). The asymptotic three-momenta of the two ejected protons 
are defined as pi and p2. We introduce relative p = Ei^Ei and center-of-mass P = pi + P2 
momenta of the ejected pair. Factorizing the y4(7,pp) and A{e,e'pp) cross sections requires 
the assumption that the sudden emission of two correlated protons only occurs when they 
reside in a relative S state llSj]. This is a reasonable approximation as investigations of the 
^^0{e,e'pp) reaction at the electron accelerators in Mainz [igI and Amsterdam 17| have 
shown that proton pairs are exclusively subject to short-range correlations when they reside 
in a relative 5* state under conditions corresponding with relatively small center-of-mass 
momenta P. In other words, the distinctive feature of the correlated two protons is that 
they are very close and moving back-to-back. With the relative S'-state assumption, the dif- 
ferential cross section is proportional to the distorted momentum distribution corresponding 



5 



5 [fill-] (5l/2-sl/2) RMSCA I 



b[fin'-I 




lfin]4 



3Rlftnl4 



FIG. 2: (color online) The function 6{R,6) for the exclusive ^^0(7, pp) cross section. In ah situa- 
tions we consider an energy transfer of 3 GeV and a three-momentum transfer q that is tuned to 
probe the maximum of the momentum distribution PniKi,n2K2 (i-^- -P=0 MeV for knockout from 
the (sl/2 - sl/2)- and (p3/2 - p3/2)-orbits and P=160 MeV for removal from the (sl/2 - p3/2)- 
orbits). We consider coplanar and symmetric kinematics, i.e. the two escaping protons have the 
same energy and polar angle 9pq, but escape from the opposite side of q. 



with the center-of-mass motion of the pair 



(p) = V 



si,S2,mi,m2 
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dR 



= dR d9S {R, 9) 
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here, = i>L+P2~^is the pair's missing momentum in the quasifree approximation. It can 
be interpreted as the center-of-mass momentum of the correlated proton pair that absorbs 



the proton. The quantum numbers (ni^i, 712^2) determine the orbits of the two correlated 
nucleons. 

TABLE I: The average effective density p probed in the various reactions. 



Reaction 


orbits 


'p{RPWIA) 


-plyRMSGA) 






(fm-3) 




{p, 2p) 


sl/2 


0.100 


0.055 


{p, 2p) 


p3/2 


0.050 


0.025 


(e, e'p) 


sl/2 


0.100 


0.086 


(e, e'p) 


p3/2 


0.050 


0.038 


(e, e'pp) 


(sl/2 - sl/2) 


0.150 


0.135 


(e, e'pp) 


(p3/2-p3/2) 


0.095 


0.075 


(e, e'pp) 


(sl/2 -p3/2) 


0.115 


0.098 



We now present the resuhs of the numerical calculations for 5{r,d) and 6{R,6). A rel- 
ativistic single-particle model was used for ^^C with parameters adjusted to describe the 
ground-state properties. The profile functions entering the «Srmsga operator have three pa- 
rameters that have been determined from the database of proton-proton and proton-neutron 



cross sections [12|]. In Fig. [T] we display the function 6{r, 6) defined in the Eq. ([T]) for proton 
knockout from the sl/2 and p3/2 orbit from a ^^C target. We compare the {p,2p) with 
the (e, e'p) result for an energy transfer of 1.5 GeV and conditions probing the maximum 
of the undisturbed momentum distribution Pnnip). The latter can be obtained by setting 
'5rmsga = 1 in Eq. ([T]). In the considered kinematics, it is clear that in the absence of nuclear 
attenuation (>Srmsga = 1 ), the upper (0° < < 90°) and lower hemisphere (90° <e < 180°) 
of the target nucleus equally contributes to 5(r, 6) and the measured signal. Moreover, the 
5{r,6) becomes equal for (e, e'p) and {p,2p). We refer to this situation as the relativistic 
plane- wave impulse approximation (RPWIA). As the ISI and FSI have the strongest impact 
at the highest nuclear densities, the RMSGA predictions for 6{r, 9) are shifted to larger 
values of r in comparison with the RPWIA ones. In addition, the contribution from the 
upper and lower hemisphere becomes asymmetric when considering attenuation. Indeed, 
the nuclear hemisphere closest to the proton detector provides the strongest contribution 
to the detected signal. The stronger the effect of attenuation the larger 5(r, 6) experiences 



a shifts in r, the larger the induced asymmetries between the upper and lower hemisphere 
and the stronger the reduction. Obviously, the asymmetry, shift and reduction occur for 
the 6{r, 9) in (e, e'p) and (p, 2p). All three effects, however, are far more pronounced for the 
{p,2p) than for the corresponding S{r,6) in (e, e'p). 

In Fig. [2] we display for the ^^C target the function 6{R, 9) defined in the Eq. ^ for 
two-proton knockout from the (j93/2 — j93/2), (sl/2 — sl/2) and (sl/2 — p3/2) orbits. Com- 
paring Figs. [1] and [2] it is clear that two-proton removal at high energies, really succeeds in 
probing the high-density regions of the target nucleus (note the different range in the radial 
coordinate r for Figs. [Hand [2]). The attenuation mechanisms induce shifts to the surface 
but the bulk of the measured strength can be clearly attributed to high-density regions in 
the target nucleus. 

In order to quantify the average densities that the various reactions can probe, we intro- 
duce js, 4] 

_ _ / drdOp if) 6 (r, 6) 
^ ~ Jdrde5 ir,e) ' 

where p (r) is the density of the target nucleus and 6 (r, 6) the function as it was defined 
in Eq. ([1]) (single-proton knockout) and Eq. ([3]) (two-proton knockout). Table [T] lists a 
systematic comparison of the computed values of p. The average density probed in the two- 
proton removal reaction from the (sl/2 — sl/2) orbits approaches the nuclear saturation 
density of po = 0.17 fm~^. We wish to stress the strong dependence on the nuclear orbit. 
Despite the strong attenuation, the {p, 2p) reaction from the sl/2 orbit can effectively probe 
higher densities than the (e, e'p) reaction from the valence p3/2 shell. For the {p, 2p) reaction 
with knockout from the sl/2 orbit the predicted effective mean density from the RMSGA 
calculations is p ~ 0.33 po- This number is almost identical to the DWIA results of Ref. ^ 
for ^^C(p, 2p) for 1 GeV incoming protons. 

In summary, we have used a relativistic framework to make a comparative and consistent 
study of the effective nuclear densities that can be probed in {p,2p), {e,e'p) and {■y,pp) 
reactions. As a representative example we have selected a carbon target and ejected proton 
kinetic energies of 1.5 GeV. We consider the results as representative for light nuclei and 
sufficiently high kinetic energies. The conclusions drawn in this work are of importance for 
ongoing and planned searches of nuclear effects at small distance scales. The (e, e'p) reaction 
has the potential to probe reasonable densities. Of all reactions considered here, the {'y,pp) 
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reaction is the one that can get closest to the deep nuclear interior. The (p, 2p) reaction is 
subject to large attenuation, but a high resolution experiment picking protons from s-orbits, 
for example, can probe densities that are of the order of 30% of po- 
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